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Abstract 

A  serpentine  flow  field  with  outlet  channels  having  modified  heights  or  lengths  was  designed  to  improve  reactant  utilization  and  liquid  water 
removal  in  proton  exchange  membrane  (PEM)  fuel  cells.  A  three-dimensional  full-cell  model  was  developed  to  analyze  the  effects  of  the  contraction 
ratios  of  height  and  length  on  the  cell  performance.  Liquid  water  formation,  that  influences  the  transport  phenomena  and  cell  performance,  was 
included  in  the  model.  The  predictions  show  that  the  reductions  of  the  outlet  channel  flow  areas  increase  the  reactant  velocities  in  these  regions, 
which  enhance  reactant  transport,  reactant  utilization  and  liquid  water  removal;  therefore,  the  cell  performance  is  improved  compared  with  the 
conventional  serpentine  flow  field.  The  predictions  also  show  that  the  cell  performance  is  improved  by  increments  in  the  length  of  the  reduced 
flow  area,  besides  greater  decrements  in  the  outlet  flow  area.  If  the  power  losses  due  to  pressure  drops  are  not  considered,  the  cell  performance 
with  the  contracted  outlet  channel  flow  areas  continues  to  improve  as  the  outlet  flow  areas  are  reduced  and  the  lengths  of  the  reduced  flow  areas 
are  increased.  When  the  pressure  losses  are  also  taken  into  account,  the  optimal  performance  is  obtained  at  a  height  contraction  ratio  of  0.4  and  a 
length  contraction  ratio  of  0.4  in  the  present  design. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Proton  exchange  membrane  fuel  cells  (PEMFCs)  are  attrac¬ 
ting  considerable  interests  as  transport,  stationary  and  portable 
power  sources.  Modeling  and  simulations  are  being  used  exten¬ 
sively  in  research  institutions  and  industries  across  the  world 
to  gain  a  better  understanding  of  the  fundamental  processes, 
so  that  the  design  cycle  can  be  shortened  significantly.  During 
the  past  decades,  many  analyses,  models  and  numerical  simu¬ 
lations  have  been  developed  [1-14]  to  study  various  transport 
phenomena  and  electrochemical  kinetics,  thus  to  gain  a  better 
understanding  and  to  develop  strategies  for  optimal  design  and 
operation  scenarios. 

Springer  et  al.  [1]  presented  a  one-dimensional  steady-state 
model  for  a  PEMFC  based  on  experimentally  determined  trans¬ 
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port  parameters.  Bernardi  and  Verbrugge  [2,3]  developed  a 
full-cell  model  to  investigate  the  mechanisms  of  water  transport. 
Fuller  and  Newman  [4]  and  Nguyen  and  White  [5]  proposed 
two-dimensional  transport  models  to  examine  the  water  mana¬ 
gement  problems  in  a  PEMFC.  Ge  and  Yi  [6]  developed  a 
two-dimensional  model  to  investigate  the  effects  of  operating 
conditions  and  membrane  thicknesses  on  the  water  transport. 
Okada  et  al.  [7]  studied  the  water  transport  at  the  anode  side  and 
gave  a  linear  transport  equation  based  on  the  water  diffusion 
and  electro-osmotic  water  drag  to  analyze  the  water  concentra¬ 
tion  profiles.  Gurau  et  al.  [8]  considered  the  variations  of  the 
concentrations  and  the  partial  pressures  in  the  gas  channels  and 
developed  a  two-dimensional  model  for  the  entire  sandwich  of  a 
PEMFC.  They  further  derived  a  half-cell  model  for  the  cathode 
side  and  obtained  rigorous  analytical  solutions  which  account 
for  the  liquid  water  content  in  the  gas  diffusion  layer  (GDL)  [9] . 
Um  et  al.  [10]  developed  a  multi-dimensional  model  to  study 
the  electrochemical  kinetics,  current  distributions,  fuel  and  oxi¬ 
dant  flows,  and  multi-component  transport  in  a  PEMFC  with 
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an  interdigitated  flow  field.  Djilali  and  Lu  [11]  focused  on  the 
modeling  of  non-isothermal  and  non-isobaric  effects  to  analyze 
the  cell  performance  and  water  transport  over  a  range  of  opera¬ 
ting  current  densities.  A  quasi-three-dimensional  model  of  water 
transport  in  PEMFCs  was  proposed  by  Kulikovsky  [12]  to  ana¬ 
lyze  the  non-linear  diffusion  of  liquid  water  in  the  membrane. 
Mazumder  and  Cole  [13,14]  developed  a  three-dimensional 
model  to  investigate  the  cell  performance  of  PEMFCs  with 
consideration  of  the  liquid  water  effects. 

The  flow  field  design  in  the  bipolar  plates  of  a  fuel  cell  is 
one  of  the  most  important  issues  for  a  PEMFC.  An  appro¬ 
priate  flow  field  design  in  the  bipolar  plates  can  improve  the 
reactant  transport,  the  thermal  and  water  management.  To  this 
end,  different  flow  field  configurations,  including  parallel,  ser¬ 
pentine,  interdigitated,  and  many  other  combined  versions, 
have  been  developed.  Many  efforts  have  been  devoted  to 
optimize  the  flow  field  design  to  improve  cell  performance 
[15-38]. 

Recently,  Soong  et  al.  [39]  proposed  a  relatively  novel  confi¬ 
guration  of  partially  blocked  reactant  channels.  They  focused 
on  the  blockage  effects  of  various  gap  ratios  and  numbers  of 
baffles  on  the  reactant  transport  and  the  pressure  losses.  They 
found  that  reducing  the  gap  size  and/or  increasing  the  num¬ 
ber  of  baffles  enhanced  the  reactant  transport.  Liu  et  al.  [40] 
developed  a  two-dimensional  model  to  examine  the  reactant 
transport  and  the  cell  performance  of  a  PEMFC  with  a  tape- 
red  parallel  flow  channel.  Numerical  predictions  show  that  the 
cell  performance  is  enhanced  by  the  tapered  reactant  chan¬ 
nel,  and  the  enhancement  is  more  evident  at  lower  voltages. 
Later,  Yan  et  al.  [41]  proposed  a  parallel  flow  channel  design 
with  tapered  heights  or  widths  to  improve  reactant  utilization 
efficiency  for  PEMFCs.  A  three-dimensional  numerical  model 
predicted  that  with  the  tapered  channel  designs,  the  flow  area 
contraction  along  the  flow  channel  leads  to  increased  reactant 
velocities;  thus,  enhancing  the  reactant  transport  through  the 
porous  layers,  reactant  utilization  and  liquid  water  removal.  The 
results  also  show  that  the  cell  performance  can  be  improved 
by  either  decreasing  the  height  taper  ratio  or  increasing  the 
width  taper  ratio.  However,  these  investigations  utilized  sim¬ 
plified  two-dimensional  models,  or  only  considered  essential 
parts  of  the  fuel  cell  in  the  computational  domain  to  reduce 
the  computational  times.  In  addition,  these  investigations  only 
considered  the  effect  of  tapered  flow  channels  on  the  parallel 
flow  field,  while  the  effect  on  a  serpentine  flow  field  has  not 
yet  been  analyzed  in  the  available  literature.  This  work  presents 
an  analysis  of  a  novel  serpentine  flow  field  design  with  redu¬ 
ced  outlet  channel  heights  or  lengths  to  improve  the  efficiency 
of  reactant  transport,  reactant  utilization,  and  cell  performance 
in  PEMFCs.  A  three-dimensional,  full-cell  model  was  develo¬ 
ped  to  analyze  the  effect  of  the  outlet  channel  height  or  length 
reductions  on  the  cell  performance  and  reactant  transport.  The 
effects  of  liquid  water  formation  on  the  reactant  transport  were 
taken  into  account  in  the  model.  The  oxygen  mass  flow  rates 
and  the  liquid  water  distributions  at  the  interface  between  the 
cathode  GDL  and  the  catalyst  layer  (CL),  and  the  local  current 
densities  are  analyzed  to  show  the  advantages  of  the  channel 
size  reductions. 


2.  Numerical  model 

Fig.  1  presents  schematics  of  the  three-dimensional  PEMFCs 
considered  in  the  present  paper.  The  cell  included  the  anode  flow 
channels,  membrane  electrode  assembly  (including  the  anode 
GDL,  anode  CL,  proton  exchange  membrane  (PEM),  cathode 
CL,  and  cathode  GDL),  and  cathode  flow  channels.  The  cell 
has  a  dimension  of  23  mm  x  23  mm  x  2.645  mm,  the  diffusion 
layer  is  0.3  mm  thick,  the  CL  is  0.005  mm  thick,  and  the  PEM 
is  0.035  mm  thick.  The  widths  and  heights  of  the  channels  and 
ribs  are  all  1  mm. 

This  paper  analyses  the  effects  of  the  reduced  outlet  chan¬ 
nel  sizes  on  the  cell  performance  and  reactant  transport.  The 
height  contraction  ratio  of  the  outlet  channel,  <$h,  is  defined  as 
8u=Hc/Hq,  and  the  length  contraction  ratio  of  the  outlet  chan¬ 
nel,  <$l,  is  defined  as  8l=Lc/Lo  (as  shown  in  Fig.  2),  where  Hc 
is  the  reduced  height  of  the  outlet  channel,  Ho  is  the  height  of 
the  inlet  channel  (1  mm),  Lc  is  the  length  of  the  reduced  part 
of  the  outlet  channel,  and  Lo  is  the  total  length  of  the  channel 
(70.1mm). 

For  convenience  of  comparison,  the  operating  conditions 
were  all  the  same  for  all  the  cells.  The  fuel  cell  temperature 
was  assumed  to  be  323  K,  the  reactants  on  the  anode  side  were 
hydrogen  and  water  vapor  with  a  relative  humidity  of  100%, 
while  the  reactants  on  the  cathode  side  were  oxygen,  nitrogen, 
and  water  vapor  with  a  relative  humidity  of  100%,  the  inlet  flow 
rate  on  the  anode  side  was  150  cm3  min-1,  and  the  inlet  flow 
rate  on  the  cathode  side  was  300  cm3  min-1 . 

A  three-dimensional  full-cell  model  was  developed  to  ana¬ 
lyze  the  electrochemical  reactions,  the  transport  phenomena  of 
the  reactants  and  the  products  in  the  cell  using  the  finite  volume 
method.  The  governing  equations  include  the  mass,  momen¬ 
tum,  species  and  electrical  potential  conservation  equations.  The 
model  assumes  that  the  system  is  three-dimensional  and  steady; 
the  inlet  reactants  are  ideal  gases;  the  system  is  isothermal;  the 
flow  is  laminar;  the  fluid  is  incompressible;  the  thermophysical 
properties  of  the  reactants  and  products  are  constant;  and  the 
porous  layers  such  as  the  diffusion  layer,  the  CL  and  the  PEM 
are  isotropic.  More  details  on  the  numerical  model  and  solution 
procedure  were  given  elsewhere  [39-42].  All  the  parameters 
used  in  the  model  are  listed  in  Table  1 . 

3.  Results  and  discussion 

The  grid  independence  was  examined  in  the  preliminary  test 
runs.  Three  non-uniformly  distributed  grid  configurations  were 
evaluated  for  a  PEMFC  with  <$h  =  0.4  and  <$l  =  0.3  at  an  operating 
voltage  of  0.3  V.  The  numbers  of  elements  in  the  X ,  Y  and  Z 
directions  were:  (I)  93  x  93  x  42,  (II)  116x  116x48  and  (III) 
139  x  139  x  54.  The  influence  of  the  elements  number  on  the 
local  current  density  is  shown  in  Fig.  3.  The  variation  of  local 
current  densities  was  5.82-12.89%  for  grids  (I)  and  (II),  and 
only  0.76-2.71%  for  grids  (II)  and  (III).  Thus,  grid  (II)  was 
chosen  as  a  tradeoff  between  accuracy  and  execution  time  in  the 
simulations. 

In  the  present  paper,  the  oxygen  mass  flow  rates,  the  liquid 
water  distributions,  and  the  local  current  densities  on  the  spe- 
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Fig.  1.  Schematics  of  PEMFCs.  (a)  Conventional  serpentine  flow  field  and  (b)  serpentine  flow  field  with  contracted  outlet  channels. 


Fig.  2.  Schematic  of  single  serpentine  flow  channel  with  contracted  outlet  chan¬ 
nel  height  and  length. 


Fig.  3.  Influence  of  the  number  of  grid  elements  on  the  local  current  densities 
at  an  operating  voltage  of  0.3  V. 
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Table  1 


Fuel  cell  parameters 


Parameter 

Value 

9.23  x  108  Am“3 

Ajtc 

1.5  x  102  Am-3 

a  a 

0.5 

1.5 

^channel 

1 

r channel 

1 

^channel 

oom2 

£gdl 

0.4 

tGDL 

1.5 

&GDL 

1.76  x  10-11  m2 

£cl 

0.4 

rcL 

1.5 

kch 

1.76  x  10-11  m2 

^Mem 

0.28 

^Mem 

Dagan  model 

^Mem 

1.8  x  10-18  m2 

cified  cross-sections  inside  the  cell  at  the  two  representative 
operating  voltages  of  0.3  V  and  0.7  V,  as  shown  in  Fig.  1,  were 
analyzed  to  understand  the  effects  of  the  various  flow  field  desi¬ 
gns  on  the  cell  performance.  The  oxygen  mass  flow  rates  and 
liquid  water  distributions  were  analyzed  at  the  interface  between 
the  cathode  GDL  and  the  CL,  with  the  local  current  densities  ana¬ 
lyzed  at  the  middle  cross-section  of  the  PEM.  The  local  transport 
phenomena  in  the  entire  cell  were  analyzed  by  plotting  all  local 
physical  quantities  along  a  specified  serpentine  channel  (Fig.  1) 
from  the  inlet  to  the  outlet. 

3.1.  Effect  of  outlet  channel  height  contraction  ratio 

The  effect  of  the  outlet  channel  height  contraction  ratio  on 
the  cell  performance  was  analyzed  at  <5h  equal  to  0.2,  0.4  0.6 
and  0.8  with  a  fixed  <$l  =  0.3.  The  conventional  serpentine  flow 
field  (without  contractions)  was  used  as  the  baseline. 

Fig.  4  shows  the  effect  of  <$h  on  the  local  oxygen  mass  flow 
rates  through  the  cathode  GDL-CL  interface  at  operating  vol¬ 
tages  of  0.3  V  and  0.7  V,  where  L  denotes  the  distance  from  the 
flow  channel  inlet.  Fig.  4  indicates  that  at  the  higher  operating 


Fig.  4.  Oxygen  mass  flow  rates  on  the  cathode  GDL-CL  interface  for  various 
($h  with  fixed  5l  =  0.3  at  operating  voltages  of  0.3  V  and  0.7  V. 


voltage  of  0.7  V,  the  oxygen  mass  flow  rates  are  almost  unchan¬ 
ged  along  the  flow  direction  and  are  almost  equal  for  the  five 
designs.  At  the  lower  operating  voltage  of  0.3  V,  the  oxygen 
mass  flow  rates  are  all  higher  than  those  at  the  higher  operating 
voltage  of  0.7  V  for  the  five  designs,  indicating  that  at  lower 
operating  voltages  the  electrochemical  reactions  are  stronger 
with  more  oxygen  consumption.  For  all  five  designs,  the  maxi¬ 
mum  oxygen  mass  flow  rates  occur  at  the  inlet  and  gradually 
decrease  along  the  flow  directions,  with  local  maximums  at  the 
turns  in  the  flow  channels  ( L-  18.5  mm,  24.5  mm,  40.5  mm  and 
50.5  mm  as  shown  in  Fig.  1).  The  inlet  oxygen  concentration 
is  the  highest,  thus  the  largest  oxygen  mass  flow  rates  through 
the  GDL-CL  interface  occur  near  the  inlet.  Then,  the  oxygen 
is  gradually  consumed  by  the  electrochemical  reaction,  so  the 
oxygen  concentration  gradually  decreases  along  the  flow  direc¬ 
tion,  and  thus,  the  oxygen  mass  flow  rate  across  the  GDL-CL 
interface  gradually  decreases.  At  the  channel  turns,  the  corners 
of  flow  channels  force  more  reactants  to  flow  into  the  GDL 
and  CL,  thus,  the  oxygen  mass  flow  rates  through  the  GDL-CL 
interface  sharply  increase  at  the  turns.  Fig.  4  also  indicates  that 
the  local  oxygen  mass  flow  rates  oscillate  in  the  transverse  flow 
channels  (L  =  1 8.5-24.5  mm  and  L  =  40.5-50.5  mm)  because  the 
anode  flow  field  channels  are  in  parallel.  Furthermore,  near  the 
flow  channel  outlets,  where  the  flow  channel  cross-sectional  area 
reduces,  as  <$h  increases  the  flow  channel  area  decreases  and 
the  reactant  flow  velocities  increase,  which  improves  the  liquid 
water  removal  efficiency,  so  the  local  oxygen  mass  flow  rates 
increase. 

Fig.  5  shows  the  effect  of  8u  on  the  liquid  water  concentra¬ 
tions  along  the  cathode  GDL-CL  interface  at  operating  voltages 
of  0.3  V  and  0.7  V.  Fig.  5  indicates  that  for  the  five  desi¬ 
gns  the  liquid  water  concentrations  at  0.3  V  operating  voltage 
are  all  higher  than  those  at  0.7  V  operating  voltage  due  to 
the  higher  electrochemical  reaction  rates.  At  the  lower  opera¬ 
ting  voltage  of  0.3  V,  the  liquid  water  concentrations  gradually 
increase  along  the  flow  directions,  with  local  minimums  at  the 
turns  in  the  flow  channels  (L=  18.5  mm,  24.5  mm,  40.5  mm 
and  50.5  mm).  In  addition,  the  liquid  water  concentrations  also 
oscillate  in  the  transverse  flow  channel  ( L-  18.5-24.5  mm  and 


Fig.  5.  Liquid  water  distributions  on  the  cathode  GDL-CL  interface  for  various 
<5h  with  fixed  <5l  =  0.3  at  operating  voltages  of  0.3  V  and  0.7  V. 
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Fig.  6.  Local  current  densities  on  the  middle  cross-section  in  the  PEM  for  various 
<5h  with  fixed  <$l  =  0.3  at  operating  voltages  of  0.3  V  and  0.7  V. 

L  =  40.5-50.5  mm)  because  the  anode  flow  channels  are  in  paral¬ 
lel.  Fig.  5  also  indicates  that  the  reduced  height  outlet  channel 
heights  efficiently  improve  the  removal  of  liquid  water  trapped 
in  the  porous  layers.  As  8u  increases,  the  reactant  flow  velocities 
increase  and  the  liquid  water  removal  efficiency  increases. 

Fig.  6  shows  the  effect  of  8u  on  the  local  current  densities 
of  the  middle  cross-section  in  the  PEM  at  operating  voltages 
of  0.3  V  and  0.7  V.  A  comparison  of  Figs.  4-6  indicates  that 
the  local  current  density  distributions  are  similar  to  the  oxygen 
mass  flow  rate  distributions,  but  are  opposite  to  the  liquid  water 
concentration  distributions.  This  is  because  lower  liquid  water 
concentrations  at  the  cathode  GDL-CL  interface  allow  more 
oxygen  to  enter  the  CL  to  participate  in  the  electrochemical 
reactions,  so  the  oxygen  flow  rates  and  the  local  current  densities 
increase. 

The  flow  field  design  affects  not  only  the  cell  performance 
but  also  the  pressure  drop  in  the  fuel  cell.  Large  pressure  drops 
in  the  fuel  cell  mean  that  more  power  is  needed  to  pump  the 
reactants.  Thus,  the  pressure  drop  is  a  significant  issue  to  be 
considered  in  choosing  the  flow  field  designs  besides  the  /-VCell 
curve.  Fig.  7  shows  the  effect  of  <$h  on  the  local  pressure  drops 
on  the  middle  cross-section  in  the  flow  channels  at  voltages  of 
0.3  V  and  0.7  V.  Fig.  7  shows  that  the  operating  voltage  has  rela- 


Fig.  7.  Pressure  drops  on  the  middle  cross-section  in  the  flow  channel  for  various 
<$h  with  a  fixed  <5l  =  0.3  at  operating  voltages  of  0.3  V  and  0.7  V. 


Fig.  8.  Polarization  curves  of  PEMFCs  for  various  <$h  with  a  fixed  <$l  =  0.3. 

tively  small  effect  on  the  pressure  drops.  At  the  lower  operating 
voltage  of  0.3  V,  the  pressure  drops  are  slightly  higher  than  at 
the  higher  operating  voltage  of  0.7  V  because  more  liquid  water 
is  produced  and  accumulated  in  the  porous  layer  and  the  flow 
channels,  which  increases  the  flow  resistances.  However,  as  <$h 
increases,  the  pressure  drops  increase  sharply  because  the  decre¬ 
ment  of  the  flow  channel  cross-sectional  areas  and  increment  of 
the  flow  resistances.  For  example,  the  total  pressure  drop  for 
<5h  —  0.6  is  about  three  times  higher  than  that  of  the  conventional 
serpentine  flow  field,  while  for  <$h  =  0.8  it  becomes  about  ele¬ 
ven  times  higher  than  that  of  the  conventional  serpentine  flow 
field. 

Fig.  8  shows  the  /-VCell  polarization  curves  for  the  various 
<$H.  Fig.  8  indicates  that  at  high  operating  voltages,  the  cell  per¬ 
formance  is  not  dependent  on  the  flow  channel  design,  while 
at  low  operating  voltages  the  cell  performance  is  significantly 
affected  by  the  flow  channel  design.  As  8u  increases,  the  limi¬ 
ting  current  density  increases  and  the  cell  performance  improves. 
Although  larger  outlet  channel  height  contraction  ratios  improve 
the  cell  performance,  the  contractions  also  significantly  increase 
the  pressure  drops;  therefore,  there  is  an  optimal  height  contrac¬ 
tion  ratio.  With  the  present  design,  the  optimal  cell  performance 
occurs  at  8u  =  0.4. 

3.2.  Effect  of  the  outlet  channel  length  contraction  ratio 

The  effect  of  the  outlet  channel  length  contraction  ratio  on 
the  cell  performance  was  analyzed  at  8 l  equal  to  0. 1 , 0.2  0.3  and 
0.4  with  a  fixed  8u  =  0.4.  Again,  the  conventional  serpentine  flow 
field  was  used  as  the  baseline. 

Fig.  9  shows  the  effect  of  <5l  on  the  local  oxygen  mass  flow 
rates  through  the  cathode  GDL-CL  interface  at  operating  vol¬ 
tages  of  0.3  V  and  0.7  V.  Fig.  9  shows  that  the  oxygen  mass 
flow  rates  at  0.3  V  is  much  larger  than  those  at  0.7  V,  indica¬ 
ting  that  the  electrochemical  reactions  are  stronger  with  more 
liquid  water  produced.  In  addition,  the  oxygen  mass  flow  rates 
along  the  flow  direction  distribute  quite  unevenly,  so  that  the 
flow  field  designs  significantly  affect  the  cell  performance  at 
low  operating  voltages.  Fig.  9  also  shows  that  at  the  lower  ope¬ 
rating  voltage  of  0.3  V,  the  distributions  of  the  oxygen  mass 
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Fig.  9.  Oxygen  mass  flow  rates  on  the  cathode  GDL-CL  interface  for  various 
<5l  with  fixed  <5h  =  0.4  at  operating  voltages  of  0.3  V  and  0.7  V. 

flow  rates  almost  coincide  at  L<  0.0426  m  for  all  the  five  desi¬ 
gns,  because  the  flow  channel  does  not  start  to  contract  until 
L  =  0.0426m.  At  L> 0.0426m,  the  oxygen  mass  flow  rates  for 
SL  =  0.4  increase  due  to  the  higher  reactant  flow  velocity.  Simi¬ 
larly,  at  L-  0.0497  m  for  <$l  =  0.3,  L  =  0.0568  m  for  <$l  =  0.2,  and 
L-  0.0639  m  for  <$l  =  0.1,  the  oxygen  mass  flow  rates  start  to 
increase.  Therefore,  the  contracted  region  of  the  outlet  channels 
increase  as  <$l  increases,  and  the  oxygen  utilization  efficiency  is 
enhanced. 

Fig.  10  shows  the  effect  of  <$l  on  the  liquid  water  concentra¬ 
tions  along  the  cathode  GDL-CL  interface  at  operating  voltages 
of  0.3  V  and  0.7  V.  Fig.  10  indicates  that  at  the  higher  opera¬ 
ting  voltage  of  0.7  V,  the  liquid  water  concentrations  for  all  the 
five  designs  are  lower  than  those  at  the  lower  operating  voltage 
of  0.3  V  due  to  the  weaker  electrochemical  reaction  rates.  The 
liquid  water  concentrations  gradually  increase  along  the  flow 
directions  and  almost  coincide  at  L  <  0.0426  m  for  the  five  desi¬ 
gns.  At  L  =  0.0426  m,  where  the  flow  channel  starts  to  contract 
for  <$l  =  0.4,  the  reactant  flow  velocity  significantly  increases, 
which  enhances  the  liquid  water  removal,  so  the  liquid  water 
concentration  for  <$l  =  0.4  is  lower  than  those  for  the  other  four 
designs.  Similarly,  at  L  =  0.0497  m  for  <$l  =  0.3,  L  =  0.0568  m  for 
<$L  =  0.2,  and  L  =  0.0639  m  for  <$l  =  0. 1 ,  the  liquid  water  concen¬ 


L  (m) 


Fig.  10.  Liquid  water  distributions  on  the  cathode  GDL-CL  interface  for  various 
<5l  with  fixed  <5h  =  0.4  at  operating  voltages  of  0.3  V  and  0.7  V. 
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Fig.  11.  Local  current  densities  on  the  middle  cross-section  in  the  PEM  for 
various  <5l  with  fixed  <5h  =  0.4  at  operating  voltages  of  0.3  V  and  0.7  V. 

trations  start  to  decrease.  Therefore,  as  <$l  increases,  which 
means  longer  flow  channel  contractions,  the  liquid  water  remo¬ 
val  efficiency  enhances. 

Fig.  1 1  shows  the  effect  of  <$l  on  the  local  current  densities  at 
the  middle  cross-section  of  the  PEM  with  operating  voltages  of 
0.3  V  and  0.7  V.  A  comparison  of  Figs.  9-1 1  again  indicates  that 
the  local  current  density  distributions  are  similar  to  the  oxygen 
mass  flow  rate  distributions,  but  are  opposite  to  the  liquid  water 
concentration  distributions.  Therefore,  increasing  <$l  improves 
cell  performance  due  to  the  higher  oxygen  mass  flow  rates  and 
the  lower  liquid  water  concentrations  in  the  porous  layers  under 
the  contracted  outlet  channels. 

Fig.  12  shows  the  effect  of  <$l  on  the  local  pressure  drops  at 
the  middle  cross-section  of  the  flow  channels  with  voltages  of 
0.3  V  and  0.7  V.  Fig.  12  indicates  that  although  the  increments 
of  <$l  promote  the  pressure  drops,  the  amount  of  pressure  drops 
are  very  small  compared  to  those  caused  by  increments  of  <$h- 
For  example,  the  pressure  drop  is  95.1  Pa  for  <$L  =  0.4,  85.4  Pa 
for  <$l  =  0.3,  75.1  Pa  for  <$l  =  0.2  and  66.6  Pa  for  <$l  =  0.1,  with 
increments  of  only  76.1%for<$L  =  0.4,58.1%for<$L  =  0.3,39.1% 
for  <$l  =  0.2  and  23.3%  for  <$l  =  0.1,  compared  to  54.0Paforthe 
conventional  serpentine  flow  field  design. 


Fig.  12.  Pressure  drops  on  the  middle  cross-section  in  the  flow  channel  for 
various  with  a  fixed  <5h  =  0.4  at  operating  voltages  of  0.3  V  and  0.7  V. 
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Fig.  13.  Polarization  curves  of  PEMFCs  for  various  <$l  with  a  fixed  <5h  =  0.4. 

Fig.  13  shows  the  I-V ceil  polarization  curves  for  the  various 
<$L,  which  indicates  that  as  <$l  increases,  the  limiting  current 
density  increases  and  the  cell  performance  improves  because 
larger  <$l  means  longer  contraction  areas  in  the  outlet  chan¬ 
nels,  which  enhances  the  oxygen  utilization  and  the  liquid 
water  removal.  Because  increments  in  the  outlet  channel  length 
contraction  ratios  only  slightly  increase  the  pressure  drop,  the 
optimal  cell  performance  occurs  at  <$l  =  0.4  with  the  present 
design. 

4.  Conclusions 

A  three-dimensional  numerical  model  was  used  to  examine 
the  effects  of  reductions  of  the  outlet  channel  flow  area  on  cell 
performance  and  local  transport  phenomena.  The  conclusions 
drawn  from  the  analyses  are: 

(1)  For  operating  voltages  greater  than  0.7  V,  the  electrochemi¬ 
cal  reaction  rates  are  lower  with  less  oxygen  consumption 
and  less  liquid  water  production,  then  the  cell  performance 
is  not  dependent  on  the  flow  channel  design.  However,  for 
lower  operating  voltages,  the  electrochemical  reaction  rates 
gradually  increase  as  the  operating  voltage  decreases.  The 
oxygen  transport  capability  and  the  liquid  water  removal 
efficiency  then  differ  for  the  various  flow  fields;  therefore, 
the  cell  performance  is  strongly  dependent  on  the  flow  field 
design. 

(2)  With  a  fixed  <$l  =  0.3,  larger  8u  causes  the  reactant 
velocities  in  the  contracted  outlet  channels  to  increase 
significantly,  which  enhances  liquid  water  removal  and 
oxygen  transport,  as  a  result,  the  cell  performance  is 
significantly  improved.  Taking  the  pressure  drop  into 
account,  the  optimal  performance  occurs  at  <$h  =  0.4  in  this 
design. 

(3)  With  a  fixed  8u  =  0.4,  larger  <$l  increases  the  outlet  channel 
contraction  region,  enhances  the  oxygen  mass  flow  rates  and 
the  liquid  water  removal  efficiency  significantly,  so  the  cell 
performance  improves. 
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